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We study the effect of gluon number fluctuations (Pomeron loops) on the impact parameter 
behavior of the scattering amplitude in the fixed coupling case. We demonstrate that the dipole- 
hadron cross-section computed from gluon number fluctuations saturates the Froissart bound and 
the growth of the radius of the black disk with rapidity is enhanced by an additional term as 
compared to the single event case. We flnd that the physical amplitude has a Gaussian impact 
parameter dependence once the gluon number fluctuations are included. This indicates that the 
fluctuations may be the microscopic origin for the Gaussian impact parameter dependence of the 
scattering amplitude. 
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I. INTRODUCTION 
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The consistent description of the impact parameter behavior of the scattering amplitude is a long standing problem. 
P^ In this vifork we discuss the influence of gluon number fluctuation on this behavior. 

In fact, some activities tovifard understanding how fluctuations change the impact parameter dependence of the 
I scattering amplitude have already started (see for example Refs. [H-Q)- Nevertheless, it is necessary to admit that we 
are still far away from a complete and consistent theory related to this subject. The most crucial difficulty related to 
^-H ■ the impact parameter dependence is the non-perturbative (soft) contribution, which should be taken into account at 
^ ' large values of the impact parameter (for reviews see [^^H and references therein) . In our approach we use a different 
J;^ technique compared to Refs. [H-I^. Namely, we calculate the rapidity dependence of the radius of the black disk in 
the fluctuation-dominated (diffusive scaling) region at high energy. 
, The gluon number fluctuations become important at very high energy. Therefore, when considering the way how 
• ' the Froissart bound may emerge based on the knowledge gathered in the small-x physics, the effects of the most recent 
, elements in the evolution and the effects of Pomeron loop, have to be taken into account. 

' In this work we focus on the consequences of fluctuations on the impact parameter dependence of the scattering 
amplitude. In Section |lT] and IIIIl we will briefly review the Froissart bound and the non-perturbative input of the 
scattering amplitude. The Froissart bound including gluon saturation effects will be studied in Section IIVI It turns 
out that the total cross section saturates the Froissart bound in the case of gluon saturation. In Section |Vl we 
will compute the impact parameter dependence of the physical amplitude including gluon number fluctuations. We 
^ ' find that the physical amplitude has a Gaussian dependence on the impact parameter, which is in agreement with 
experimental measurements. We also calculate the radius of the black disk including gluon number fluctuations and 
find a unique rapidity dependence coming from fluctuations. Further, in Section IVll we calculate the slope parameter 
B. The summary is given in Section [VIII 



II. UNITARITY AND FROISSART BOUND 



In high energy scattering processes, the total cross section for the reaction 12 — > n particles is 

Here |pi| is the magnitude of the initial center-of-mass frame three momentum. It is well known that the probability 
conservation in the scattering processes requires the scattering S matrix to be a unitary matrix SS'' = 1. Unitarity of 
S matrix provides a simple way to derive total cross sections from the forward(0s — 0) elastic scattering amplitude. 
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which is known as the optical theorem. For any orthonormal states \i) and {j\, one has 



S,.^{J\SS^\^)^J2(3\S\f){f\S^\^) (2) 
/ 



where we have used the completeness relation 



We define the T matrix as 5* = 1 — iT, then the miitarity condition requires that 

- 0-|rt|z) = i27r)Hj25\pf - P')OlTt|/)(/|T|z). (4) 
/ 

In the case of j = i, where the final state is the same as the initial state, wc obtain the optical theorem 

^ Im(z|rK) (5) 



~ 2\p,\V-s 
with 

2lm(z|r|z) = ImA{s,t^ 0) = J2{2n)SS\pf - P^)|(/|r|z)p, (6) 

/ 

where A{s,t) is the elastic scattering amplitude, s and t are the center-of-mass energy squared and the momentum 
transfer squared, respectively. 

In high energy physics, the two particles scattering amplitude A(s,cos6's) can be expanded in the partial-wave 
series, 

oo 

^(s,cos6',) = 1677^(2; + 1) A; (cos 61s), (7) 

1=0 

where P/ (cos 0s) is the Legendre polynomial of the first kind, and 9^ is the s channel scattering angle in the center- 
of-mass frame. The partial- wave amplitude Ai(s) can be written in terms of a real phase shift Sg and an inelastic 
threshold rji 

2^^) ' 

where p(s) = 2|pi|/y^ with our choice of normalization, and unitarity requires that < 77/ < 1. Ai{s) will be 
exponentially small for 

; > aA^-l^/iln(s) (9) 

and the scattering amplitude (O may be truncated at this value 7] . With the unitarity constrain < r^; < 1 and ([8|) , 
we can get 



\Ms)\ 



2ip{s) 



where p{s) — > 1 as s — > 00. We know that the Legendre polynomial of the first kind |P;(cos6's = 1)| < 1. So, for 
large s 

|A(s,cos6's = 1)1 < ^(2/ + l) (11) 

1=0 

with [ij 

= aA/"iyiln(s). (12) 
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Performing the summation over I in ([TT|) gives 

\A{s,cos9s = 1)1 < constat x sln^(s). (13) 
Substituting it into (O, the Froissart bound is 

CT*°* < constant x In^(s) (14) 

where we have used 

IPiI's = (Pi • P2)' - mlml = l[s - (mi + m.2f][s - (mi - ma)']. (15) 

Here, Pi and P2 are the magnitude of the initial center-of-mass frame four momentum in the two-body scattering 
process 1 + 2 — > 3 + 4, the toi and m2 are the mass of particle 1 and particle 2. The Froissart bound is one of the 
outstanding results of the analytic S matrix theory. This bound has been derived in 1961 by Froissart Q assuming that 
the two particles scattering amplitude has uniformly bounded partial wave amplitudes and satisfies the Mandelstam 
representation with a finite number of subtractions. The Froissart bound expresses that the hadronic total cross 
section can not rise faster than constant x In' s. 



III. THE NON-PERTURBATIVE INPUT 



In many practical applications, it is too complicated to perform calculations of the scattering amplitude keeping 
precisely the information about the impact parameter dependence, since it is related to the non-perturbative physics. 
In order to simplify the situation, one considers the scattering at fixed impact parameter and then introduces the 
knowledge about the impact parameter dependence through some profile function, which we will denote by S{b). 
Usually, the following two ansatze are used as a non-perturbative input: 

1. The scattering amplitude expressed as the product of the scattering amplitude at fixed impact parameter times 
the profile function S{b) 



T(Y,r,b) = T{Y,r)-S{b). (16) 
Such factorization form is usually used in the region of large values of the impact parameter b. 

2. The second one is mostly inspired by the numerical study of BK equation with a modified BK kernel in which 
the kernel of the BK integral equation is regulated to cut off infrared singularities [l3|. The impact parameter 
dependence is introduced through the saturation scale, Qs{Y,b) — Qs{y) ■ S(b), and consequently for the 
scattering amplitude we have: 



T{Y,Q,b) = T{Q,Qs{Y,b)) = T{Q, Qs{Y,b ^ 0) ■ S{b)) . (17) 

In both cases, the impact parameter profile function typically has the exponential behavior S{b) = e"^™"'' at large 
distances b ^ Rq, where Rq is the typical radial size of the hadron under consideration and Rq increases as A^^^ for 
a nucleus with atomic number A. We use such an exponential fall-off at large impact parameter as a non-perturbative 
initial condition at low energy. 

IV. SINGLE EVENT AMPLITUDE 

In the geometric scaling region and in the fixed coupling case, the dipole-hadron scattering amplitude reads 

T{Y,r,b)c,{r'Ql{Y)y^^ -Sib), (18) 

where the saturation momentum is 

Ql{Y) = g^e"'^^^ (19) 
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with the arbitrary reference scale Qo (Qo ~ C'(Aqcd)) and with the (6) giving the impact parameter dependence. 
Note that here factorization is assumed, which is the case as in Refs.|10l. Eq. (IT51) shows geometric scahng with 
the anomalous dimension 



7 = 1 - 7, ~ 0.37. (20) 

Sib) ~ (21) 

TiY,r,R)^K^O{l) (22) 

si; h'' -'■(!))- 



Now, with the non-perturbative input 
at large 6, one obtains from the condition 
the "black disc radius" 



Eq. (p3)) gives the standard result given in the literature . We have gone through such a detailed derivation of (|23)) 
since one of the main purposes of the present work is to show how Eq. (j23p is modified once gluon number fiuctuation 
effects are included. 

The resulting dipole-hadron cross section saturates the Froissart bound 



cr*°* ^ 2 J (fbT{Y,r,b) 



2ttR^ 



2^7.' ( - ^ , fQ 
^[-^'^sY-ln 



7r7^ / ujus \ . 2 



2 V m 



In s (24) 



with Y = ln(s/Q2). 

We would like to note that there was some controversy between [ll| and [H, [l^- The authors of Refs.p^.[l3j claim 
that the exponential fall-off with b of the initial distribution should replace by a power law fall-off due to perturbative 
nature of BK equation, which would be then too slow to satisfy the Froissart bound. While Ref . [lH points out that 
the BK equation will preserve the exponential tails at very high energy due to a quasi-locality of the BK equation. 
The numerical study of the BK equation with a modified BK kernel shows that the exponential tail is preserved when 
the kernel of the BK equation is properly regulated in the infrared fl(J]. 

It is easy to check that both ansatze in Eq. ([TE)) and Eq. pT)) for the single event amplitude lead to the similar 
result in Eq. 

V. INCLUDING GLUON NUMBER FLUCTUATIONS 

After including fluctuations one has to distinguish between the even-by-event amplitude and the averaged (physical) 
amplitude. They can be explained by considering the evolution of a hadron from y = wp to y = Y which is probed 
by a dipole of size r, giving the amplitude T{r,Y). The evolution of the hadron is stochastic and leads to random 
gluon number realizations inside the hadron at Y, corresponding to different events in an experiment. The physical 
amplitude, T{r, Y)^ is then given by averaging over all possible gluon number realizations/events, T{r, Y) = {T{r, Y)) 
, where T{r,Y) is the amplitude for the dipole r scattering off a particular realization of the evolved hadron at Y . 

Based on the high energy QCD/statistical physics correspondence, we can write [T^ 

r(p,p,(r,6)) - {T{p,p,{Y,b))) = / dps{Y,b)T{p^ p,{Y,b))P{Ps{Y,b)~ {p,{YM) (25) 

where we have used p = \n{Q'^/QQ) and ps{Y,b) = \n{Ql{Y,b) / Qq) and we have assumed the dependence of the 
scattering amplitude on the impact parameter through saturation scale [l5|. The probability distribution of /9s(y, 6) 
is argued to have Gaussian formfl6|. 



P(,.(y,5)).-^exp 



(p.(y,&)-(p,(y,6)))^ - 

DY ^ ' 
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and the single scattering amplitude T{p ~ PsiX: b)) is 



rrf (VKW ioT p<pg{Y,b) 

T{p,p.s[Y,b)) ^ < ,^ , ^ (27 

exp [-js{p - Ps{Y, b))] for p > ps{Y, b). 



It is easy to show that in the diffusive scaling region, <C p — {ps{Y, b)) ^ 7scr^, 

p-{Ps{Y,b)) ^'^^ 
where we have used = DY with D being the diffusion coefficient. Now, with the non-perturbative input 

{p,{Y, b)) ~ p,(r, b) ~ uasY - 2m,fe, (29) 
one can easily see that the exponential decrease with b in the single event case is turned into a Gaussian b dependence 

Ijy 



{T(Y,r,b)) ~ : ■ e-(p-""^^+2'"''')V^y 

p~ {ps{Y,b)) 

I 4'"Jf>^ 

CX -—e- — . (30) 

This consequence of fluctuations seems to be supported by the experimental observations [l7l - [l9j since, say for pp 
collision, 

^ ^ e-^l*l (31) 
at 

which after a Fourier transform gives 

S{b)^e-&. (32) 

Second consequence of fluctuations is that the factorization is broken, see Eq. (pS)) . 

(T(y,r,6)) ^/(y,r).5(&) (33) 

as compared to the single event amplitude in Eq. (fT8)) where 

(T(y,r,&)) =/(y,r).5(6). (34) 

Third consequence of fluctuations is that also the averaged amplitude in Eq. (1501) satisfies the Froissart bound. 
Namely, from the condition (T) = k ~ 0(1) (but k < 1), 

n ^ ^ ■ e-iP-^'^.y+^m.nf/DY (35) 
p- {ps{Y,R)) 

which after taking the logarithm on both sides reads (k' close to one) 

(p-(p,(y,i?)))2 , / Vdy 



DY ''""[p-ip^Y^R))) 
and is fullfilled if 

P-{p,{Y,R))^cVdY (37) 



with the coefficient c of order v k' . 

Now, with {ps{Y, R)) ~ wusY — 2m7ri?, one obtains 



p - ujctsY + 2m.^R = cVDY 

2m^ V ' - - VOo 



i? ^ ( wa,r + c^/~DY - In ( ^ ) ) . (38) 



As compared to Eq. (|23l) . this equation taking fluctuations into account contains a new term ^/iJY . 

So, including fluctuations and the impact parameter dependence in the way presented here seem to lead to reasonable 
results. However, the whole discussion is valid only in the flxed coupling case. We have already calculated the dipole- 
hadron scattering amplitude by considering both gluon number fluctuations and runing coupling effect [20[; however 
we do not know how to take into account the impact parameter dependence of the scattering on top of gluon number 
fluctuations and runing coupling effect. This will be considered in our next work. 
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VI. PHENOMENOLOGICAL APPLICATIONS AND ESTIMATION OF THE SLOPE PARAMETER B 



It is well known from numerous hadronic scattering experiments that 

cx e-^l*l. (39) 

at 

where t is the squared four momentum transfer between the projectile and target. The i— slope B tends to a universal 
value determined by the proton shape alone {2l|. From the experimental measurement of the distribution of the 
vector mesons, the effective slope B is found to be i? = AGeV~^ p7l - [l9| . 
In order to study the i— slope B, we take the Fourier transform of Eg . ([551) : 

Se.pib) (X e-^ . (40) 

Now from comparison of factors in the exponent of Eq.(|30|) with Eq. (|40p . we can immediately see that 

B = 



8ml 
DY 



where D is the diffusion coefficient. The value oi D ^ 0.325 is determined by fitting the HERA data with the color 
glass condensate model plus the gluon number fluctuations '55, . Note that the B increases logarithmically with 
the center of mass energy s, which is in agreement with the Rcgge theory. This is the phenomenon known as the 
shrinkage of the diffraction peek in Regge theory, which can be interpreted as an increase of the interaction radius 
Rint ^ Vlns. With the reasonable values TOtt = 0.14GeT^, ^/s — 90GeV and = lOGeV'^, one obtains values 

D / s 

B = ^ In 

8ml \Q 

~ 13.8GeV-^. (42) 



VII. SUMMARY 



The main results of this work can be summarized as follows: We have argued that the impact parameter behavior 
of the scattering amplitude in the presence of fluctuations has Gaussian-like behavior. Such behavior is in agreement 
with various phenomenological models. This indicates that fluctuations may be the microscopic origin for the Gaussian 
behavior. 

Further, we have shown that the factorization of the impact parameter of the scattering amplitude is lost once the 
gluon number fluctuations are included. 

We calculated the rapidity dependence of the radius of the black disk in the fluctuation-dominated (diffusive scaling) 
region at high energy. We found that, due to fluctuations, the growth of the radius of the black disk is enhanced by 
an additional (proportional to square root of rapidity) term. 
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